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ABSTRACT 

To better understand the physical properties of accretion disks in high-mass star formation, we 
present a study of a dozen high-mass accretion disk candidates observed at high spatial resolution 
with the Australia Telescope Compact Array (ATCA) in the high-excitation (4,4) and (5,5) lines of 
NH3. All of our originally selected sources were detected in both NH3 transitions, directly associated 
with CH3OH Class II maser emission and implying that high-excitation NH3 lines are good tracers of 
the dense gas components in hot-core type targets. Only the one source that did not satisfy the initial 
selection criteria remained undetected. From the eleven mapped sources, six show clear signatures 
of rotation and/or infall motions. These signatures vary from velocity gradients perpendicular to the 
outflows, to infall signatures in absorption against ultracompact Hii regions, to more spherical infall 
signatures in emission. Although our spatial resolution is '--^lOOOAU, we do not find clear Keplerian 
signatures in any of the sources. Furthermore, we also do not find flattened structures. In contrast 
to this, in several of the sources with rotational signatures, the spatial structure is approximately 
spherical with sizes exceeding 10* AU, showing considerable clumpy sub-structure at even smaller 
scales. This implies that on average typical Keplerian accretion disks - if they exist as expected - 
should be confined to regions usually smaller than 1000 AU. It is likely that these disks are fed by the 
larger-scale rotating envelope structure we observe here. Furthermore, we do detect 1.25 cm continuum 
emission in most fields of view. While in some cases weak cm continuum emission is associated with 
our targets, more typically larger-scale Hii regions are seen offset more than 10" from our sources. 
While these Hii regions are unlikely to be directly related to the target regions, this spatial association 
nevertheless additionally stresses that high-mass star formation rarely proceeds in an isolated fashion 
but in a clustered mode. 

Subject headings: ISM: kinematics and dynamics - Stars: rotation - Stars: formation - Stars: early- 
type Stars: individual: G305.21+0.21, G316.81-0.06, G323.74-0.26, G327.3-0.6, 
G328.81-h0.63, G331. 28-0.19, G336. 02-0.83, G345.00-0.22, G351. 77-0.54, G0.55- 
0.85, G19.47-0.17, IRAS18151-1208 - Techniques: interferometric 



1. INTRODUCTION 

The characterization of accretion disks around 
young high-mass protostars is one of the main un- 
solved questions in h igh-ma ss star formati o n re- 
search (iBeuther et all l2007al: ICesaroni et"al] 120071 : 
IZinnecker &: Yorkell2007T ) . The controversy arises around 
the difficulty to accumulate mass onto a high-mass 
protostar when it gets larger than 8M0 because the 
radiation pressure of the growing protostar m ay be 
stron g enough to revert the g as inflow (e.g., iKahnI 
[1971 IWolfire fc Cassinellil [198711 . Different ways to 
circumvent this problem are proposed, the main are 
disk accretion from a tur b ulent gas and dust core 
(e.g., [Jiiin a & Adam^ [mi lYorke fc Sonnhalteil [200l : 
iMcKee fc Tan 200;|), competitive accretion and poten- 
tial (proto) stellar mergers at the de nse centers of evolv- 
ing h i gh-mass (proto)clusters ( e.g., iBonnell et al.|[200l 
120071 : iBallv fc Zinneckeil 120051 ). and ionized accretion 
fiows continuing throu gh th e hypercompact Hii region 
phase (e.g., Keto 2003, 20071 ). 

Over recent years, much indirect evidence has been 
accumulated that high-mass accretion disks do exist. 
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The main argument stems from high-mass molecular 
outflow observations that identify coUimated and en- 
ergetic outflows from high-mass protostars, resembling 
the properties of known low-mass star formation sites 
(e.g., [B euthcr eFall l2002at IBeuther fc ShepherdI l200l 
lArce et al.ii2007f ). Such coUimated jet-like outflow struc- 
tures are only explainable if one assumes an underly- 
ing high-mass accretion disk that drives these outflows 
via magneto-centrifugal acceleration. From a theoret- 
ical modeling approach, recent 2D and 3D magneto- 
hydrodynamical simulations of high-mass collapsing gas 
cores result in the format ion of high-mas s accretion 
disks as well (lYorke fc S onnhalterl l2002l : IKrumholz et al.l 
I2007bl: [Kratter fc Matzner .20063. Although alternative 
formation scenarios are proposed, there appears to be a 
growing concensus in the high-mass star formation com- 
munity that accretion disks should also exist in high-mass 
star formation. However, it is still poorly known whether 
such high-mass disks are similar to their low-mass coun- 
terparts, hence dominated by the central protostar and 
in Keplerian rotation, or whether they are perhaps self- 
gravitating non-Keplerian entities. 

While the indirect evidence for high-mass accretion 
disks is steadily increasing, direct observational studies 
are largely missing. This lack of observational evidence 
can be attributed to two main reasons. The first is the 
clustered mode of high-mass star formation and the typ- 
ically large distances, hence spatially resolving and dis- 
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entangling such structures is a difficult task. The sec- 
ond difhculty is to choose the right spectral line tracer 
which allows unambiguous identification and characteri- 
zation of the disk structure. Many spectral lines are ei- 
ther optically thick (e.g., CO, HCO+, CS), chemically 
difficult to interpret (e.g., HC3N, HNCO) or excited 
in the envelope and disk which causes confusion (e.g., 
HCN, CH3CN). To overcome these problems, we used the 
Australia Telescope Compact Array (ATCA) at 1.2 cm 
wavelengths including the most extended baselines (re- 
sulting in a spatial resolution < 1"), and we aimed at 
the highly excited NH3(J,K) inversion lines (4,4) and 
(5,5). NH.3 is kn own to be a dense core tracer (e.g., 
IZhang et al.l Il998l ). and the high (J,K) inversion lines 
with excitation temperatures (E'lowcr) of 200 and 295 K, 
respectively, should only be excited in the innermost and 
warm region close to the central protostars. Similarly, 
iQsorio et all (|2009D also modeled the NH3(4,4) emission 
of the collapsing hot core G31.41. Furthermore, radi- 
ation transfer calculations for 3D hydro-simulations re- 
vealed that the 1.2 cm band regime should be particu- 
larly well suited for such studies because the inner disk 
regions m ay be optically thick at frequencies above about 
100 GHz (|Krumholz et al.ll2007a[ ). This may make future 
high-mass disk studies of the innermost regions at shorter 
wavelengths with ALMA difficult. Therefore, studying 
these lines at high angular resolution with the ATCA 
allows us to penetrate deeply into the natal cores and 
study the physical properties of the predicted high-mass 
accretion disks. 

Over the last few years there have been a lot of "trial 
and error" approaches for high-mass disk studies, but 
no consistent investigation of a larger sample is pub- 
lic so far. Since the above outlined approach has been 
proven very successful in the recent ATCA high-(J,K) 
NH3 study of IRAS 18089-1732 (BcutheriLWalsh 2008), 
we are now aiming at a source sample of twelve promis- 
ing disk candidates mainly identified by previous lower 
resolution NH3 studies of 41 sources using the ATCA 
([Longmore et al.l l2007af ) and 60 sources using Mopra 
(Walsh et al. priv. comm.). The proposed sources com- 
prise the best high-mass-disk-candidate sample for the 
southern hemisphere to date (Sec. [3|). 

2. OBSERVATIONS 

Our sample of twelve sources (see Table [1} was ob- 
served during six consecutive nights between July 1st and 
7th 2008. We always observed two sources per night in a 
track-sharing mode cycling between the gain calibrators 
and sources. Table [T] lists the corresponding gain cali- 
brators (phase and amplitude) for each pair of sources. 
Simultaneously, we observed the NH3(4,4) and (5,5) in- 
version lines with the frequencies of the main hyperfine 
components at 24.139 and 24.533 GHz, respectively. The 
phase reference centers and velocities relative to the lo- 
cal standard of rest (wisr) are given in Tabled] Bandpass 
and flux were calibrated with observations of 1253-055 
and 1934-638. The spectral resolution of the observa- 
tions was 62 kHz corresponding to a velocity resolution 
of ~ 0.8 km s^^. The observations were conducted in the 
1.5D configuration including antenna 6 which resulted 
in a maximum baseline length of 4.3 km. Depending on 
the source structure (e.g., compact versus extended) and 
strength we applied different weightings (robust values 



between -2 and 2) for each source and line. To better 
trace faint and/or extended features, we occasionally ex- 
cluded antenna 6 - and hence the longest baselines - from 
the imaging process (Table [2]). The synthesized beams 
and rms are given in Table [2l 

3. SAMPLE 

The source sample was largely identified by previous 
lower resolution NH3 studies of 41 sources using the 
ATCA (jLongmore et al.ll2007aD and 60 sources using Mo- 
pra (Walsh et al. priv. comm.). All sources were se- 
lected based on strong NH3 (4,4) and (5,5) emission, they 
show outflow signatures and are prominent in other dense 
core tracers (e.g., Purcelf et al. 2006, 2009). IRAS 18151- 
1208 does not follow these identification criteria but was 
selected because of recent disk-like structures observed 
in the 1.3 mm continuum emission (Fallscheer et al. in 
prep.). The outflow orientations, which have to be per- 
pendicular to the expected disks, are known for a major- 
ity of the sources. The studied sources comprise the best 
high-mass-disk-candidate sample for the southern hemi- 
sphere to date. The sample size (12 sources) was chosen 
because it doubles the number of e xisting disk candidates 
as listed by Cesar onl et al.l ()2007f ). which were observed 
heterogenously by different groups, with different tracers 
and different selection criteria. In contrast to that, these 
new data provide a homogeneous dataset which is easier 
to interpret. Investigating such a large sample is the only 
way to characterize high-mass disk properties in a gen- 
eral way, important for a comprehensive understanding 
of high-mass star formation. 

4. RESULTS 

We will first present the observational results for each 
source individually and then put them into a general con- 
text in section [S] 

4.1. Results for individual sources 
4.1.1. G305.21+0.21 (IRAS 13079-6218) 

This region exhibits linear maser featu res with an 
appro ximate northeast-southwest direction (iNorris et al. 
Il993fl similar to the H2 emission features by iDe Buizer 
(|2003f l. We clearly detected the NH3(4,4) and (5,5) fines 
including their hyperfine structure components, and as 
shown in Figure [U the corresponding intensity- weighted 
velocity maps (1st moment maps) exhibit a velocity 
gradient in the northwest-southeast direction, approxi- 
mately perpendicular to the maser and H2 signatures. 
The extent of the velocity structure is approximately 2" 
corresponding at the adopted distance of 3.5 kpc to a size 
of '^7000AU. Are these NH3 signatures due to rotation 
from an infalling/rotating envelope and/or an embed- 
ded high-mass accretion disk? The pv-diagrams of the 
NH3(4,4) and (5,5) lines in Figure[2]also clearly show the 
velocity gradient, however, a profile similar to a Keple- 
rian disk is hardly discernable. Hence, it is more likely 
that the NH3 structure corresponds to a large rotating 
and infalling envelope that may feed a real Keplerian ac- 
cretion disk at its center. 

Figure[T]also presents the intensity-weighted line width 
maps (2nd moments) . The line width distribution of the 
(4,4) transition exhibits several positions of increased line 
width. This is likely due to a clumpy sub-structure of the 
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TABLE 1 

Observational parameters I 



Source 


R.A. 


Dec. 


^Isr 


d!> 


cal 


comment 


Refs. 




J200.0 


J2000.0 


km 


kpc 








G305.21+0.21 


13:11:13.77 


-62:34:41.2 


-38.3 


3.5 


1352-63 


Im NE-SW, H2 NE-SW 


1,2,15 


G316.81-0.06 


14:45:26.90 


-59:49:16.3 


-38.7 


2.7 




Im N-S, No H2,GF NNW-SSE, 


2,4,5,13 














7MM NNW-SSE? 




G323.74-0.26 


15:31:45.80 


-56:30:49.9 


-49.6 


3.3 


1613-586 


Im SW-NE or E-W?, H2 E-W?, GF ? 


1,4,5,6 


G327.3-0.60 


15:53:09.29 


-54:36:57.5 


-46.0" 


3.1/11.2 






7 


G328.81+0.63 


15:55:48.44 


-52:43:06.0 


-41.5 


3.0 


1613-586 


Im E-W?, H2 E-W, SiO E-W 


1,2,3,4 


G331. 28-0.19 


16:11:26.90 


-51:41:56.6 


-88.1 


5.4 




Im NNW-SSE or WNW-ESE?, H2 E-W 


1,2,3 














GF NE-SW, SiO NNE-SSW, 7mm E-W? 


4,5,13 


G336. 02-0.83 


16:35:09.30 


-48:46:47.0 


-48.5 


3.6/12.0 


1600-44 


Im N-S 


4 


G345. 00-0.22 


17:05:10.90 


-41:29:06.6 


-26.8" 


2.9/13.5 




Im E-W? 


4 


G351. 77-0.54 


17:26:42.57 


-36:09:17.6 


1.2 


2.2 


1714-336 


Im NE-SW or N-S?, CO NE-SW 


1,4,8 














OH N-S, H2O ring or NE-SW? 


16,17,18 


GO.55-0.85 


17:50:14.53 


-28:54:30.7 


18.0" 


7.7/9.4 






4 


G19.47-0.17 


18:25:54.70 


-11:52:34.1 


19.7 


1.9 


1829-106 


CO NNE-SSW, GF ? 


5,13 


118151-1208 


18:17:58.24 


-12:07:24.5 


32.8 


3.0 




H2/CO NW-SE, dust NE-SW, CH3OH 


8,9,10,11,14 



" flsr taken from Mopra spectra (Walsh et al. priv. comm.). 



* Most distances were taken from the hterature. For those where we did not find distance references, we calculated the kinematic near and 
far distances using the Galactic rotation curve bv lBrand fc Blitzj 1 IT9M) . 

Comments: Im — > linear maser with orientation, H2/SiO/CO/7mm outflows with potential orientation, GF — > "Green fuzzy" Spitzer 4.5 /xm 
elongation, dust continuum orientation, CH3OH for the last source just states Class II CH3OH maser detec tion. "?" denotes uncert ainty. 
(1) Norrisctal. (1993), (2) Dc Buizcr (2003), (3) Dc Buizcr ct al. (2008), (4) Walsh ct al. (199|), (51 ILongmore et al.l (l2007alV (6) 
[Walsh et al 120021 (7) Caswell ct al. (,1995 ). (8) ieurini ct al. (2008), (9) Sridharan ct al. (200l), (101 IBeuther et al.l (I2002bl). Clll 
IPavis et alT (120041'). (12) Fall s cheer e t al. (in prep.), (13) Longmorc ct al. (2009), (14) Longmorc & Burton (2009), (15) IBeuther et al.l 
Il2002d) . (lel lLongmore et al.l (I2007bl) . (ITl lFish et all 1200 5). (18) Forstcr ct al. (1990), (19) Zapata ct al. (2008), (20) Longmore et al. (in 
prep.) 

core and the very high optical depths of the hne (see also 
Fig. [30] for example spectra). The line width distribution 
of the (5,5) transition is simpler with a clear line width in- 
crease toward the center. This indicates that the higher 
excited line (ii^iower(NH3(5, 5)) = 295 K) probes deeper 
into the core center better allocating the position of the 
central protostar. Such a central line width increase is 
also necessary in the picture of a rotating, infalling en- 
velope with a central disks. 

We also detect an Hii region (Fig. [3]), however, that is 
located ~ 50" to the east and unlikely to be associated 
with our region of interest. 

4.1.2. G316.81-0.06 (IRAS 14416-5937) 

This region exhibits Class II CH3QH maser s aligned 
in a north-south direction ([Walsh et al.l llQQST ) as well 
as elongated Spitzer 4.5 /xm "green fuzzy", extended Ks- 
band emission (likely from shocked H2) and 7 mm contin- 
uum emission with an orien tation approximately NNW- 
SSE ( Longmore etTaLl 120091 : iLongmore fc Burto'nl [2001 ) . 
While the Spitzer 4.5 ^m emission is usually attributed to 
shocked H2 from an underly ing outflow/jet system (e.g., 
iNoriega-Crespo et al.l |2004[ ) , 7 mm continuum emission 
could in general be attributed to jet or disk emission 
(e.g.. lZapata et al.l[2005l : lArava et al.ll2007t ). 

For this region, we do detect the NH3(4,4) and (5,5) 
lines less strongly than, e.g., toward the previously dis- 
cussed region G305. 21-1-0. 21, and we had to image them 
with a lower spatial resolution (see Table [21). However, 
the structure is clearly elongated approximately in a 
north-south direction (Fig. [4]). While the Spitzer 4.5 yum 
emission co-alignment with the NH3 features would sug- 
gest an association with the outflow, this appears less 
likely for spectral lines with high excitation temperatures 
(exceeding 200 K, see gl]). It appears that for G316.81- 
0.06 our NH3 data are inconclusive whether they can be 
associated with a rotating structure or not. Similarly, 



with the given low signal-to-noise ratio the 2nd moment 
line width distributions are also inconclusive with regard 
to the location of the line width maximum. 

We do also detect extended 1.2 cm continuum emis- 
sion from an Hii region, however, that is located approx- 
imately 25" to the west and hence not directly associated 
with our target source (Fig. [5]). 

4.1.3. G323.74-0.26 (IRAS 15278-5620) 

The Class II CH3OH maser features do not follow 
a clear trend, but most emission peaks a ppear to be 
align e d mainly in an ea st-west direction ( Norris et al. 
Il993t IWalsh et all Il998l) . Furthermore, IWalsh et al 
(|2002D identify a H2 outflow structure in approximately 
east-west direction. Spitzer exhibits green fuzzy emission 
in the 4.5 /im band as well, but it is difficult to associate 
an obvious outflow direction with that. 

Our high-excitation NH3 observations also show no ob- 
vious trends, and the emission features are even slightly 
different in the (4,4) and (5,5) hne (Fig. [g. While the 
(4,4) emission is elongated in approximately southeast- 
northwest direction, it does not show an obvious veloc- 
ity gradient. In contrast to that, the (5,5) emission is 
more compact but shows no clear velocity structure ei- 
ther. Since the signatures do not coincide in both lines 
and the outflow identification is not unambiguous either, 
we refrain from further interpretation. 

The 2nd moment line width distributions exhibit the 
line width maximum close to the main group of Class II 
CH3OH maser features, hence indicating that the center 
of gas infall and star formation activity, and hence the 
location of the main protostar is likely close to that. 

The 1.25 cm continuum shows a tentative 5(t = 
1 mJy beam"^ emission feature -1-11.1"/ — 7.5" offset from 
the phase center. 

4.1.4. G327. 3-0.60 
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TABLE 2 

Observational parameters II 







Rf>t1 TYl 


IfT rms'^ 

mJy 


Peal^ 

niJy 


CjoUo.zi+U.zi 


M U f A A\ 

IN 113(4,4 ) 


U.ol X U.ob 


2.5 


29 


i^Qnc; 01 1 n 01 
CjoUo.zi+U.zi 


iN 113(0,0} 


U.oU X U.OO 


2.4 


31 


CjoUo.zi+U.zi 


cont. 


/I 
z.z X Z.4 


0.8 




Golb.oi-U.Ub 


JN 113(4,4) 


"1 00 1 OCT 

1.00 X l.zo 


3.7 


36 


Qi ft CI n rift 


IN II3 (0,Oj 


ni \y ^ /II 
Z.Ul X 1.41 


ft 
Z.D 


QQ 
00 


Qi ft CI n rift 


cont. 


1 ft/1 \y 1 ni 
1.D4 X l.Ul 






Go2o.74-U.zb 


MTU (a A\ 

JN 113(4, 4) 


n Oil; > , r\ cr /I 

U.ob X U.o4 


2.0 


19 


r^QOQ 7/1 n Oft 


IN II3 (0,0 J 


U.ol X U.Ol 


1 
z. 1 


Oft 

ZD 


r^QOQ 7/1 n Oft 


cont. 


fti \/ 1 f^fi 
Z.ol X l.OD 


n 
U.Z 




Goz7.o-U.dU 


M TU (a A\ 

JN 113(4,4) 


n nn > , n cr o 

U.9U X U.OO 


2.6 


49 


r^Qov Q n ftfit' 


INII3 (4,4 ) 


nc \/ 1 1 1 
Z.Uo X 1.11 


1 
Z. 1 


IZvS 


QOV Q n ftfl 


IN II3 (0,0) 


n Qo \/ n t;o 
u.oy X U.OZ 


Q 

Z.o 


KA 

04 


r^Qov Q n ftfifa 


1N113 (0,0) 


m \y 1 nc 
Z.Ul X l.Uo 


ft 
Z.D 


1 /I n 
14U 


(^Qov Q n ftfl 


cont 


ft/1 \/ 1 Al 
Z .04 X 1.4/ 


1 c 
1.0 






M XJ (A A\ 

IN 113(4,4 ) 


n fto \x n /I /I 
U.DO X U.44 


2.0 


— 

—82 


r^QOfi 01 _Ln ft"? 


IN rl3 (0,0) 


n ft7 \/ n /I 
U.D / X U.40 


1 c 
1.0 


— 01 


Vjozo.oi-rU.Oo 


cont. 


n fti V n zll 


z.o 




LjOo1.ZS-U.1o 


MXJ (a A\ 

IN 113(4,4 ) 


U.yl X U.OO 


1.5 


22 


1_jOvS1.Zo-U.1o 


IN II3 (O,0) 


n fin \y n c^v 
U.yU X U.O/ 


1 c; 
1.0 


1 ft 
ID 


LjOo1.Zo-U.1o 


cont. 


n ft Q n /I 
U.Do X U.4 


0.3 




GcioD.Ui-U.OJ 


■\THJ (A A\ 

IN 113(4,4) 


on > , 1 00 

Z.2U X l.zo 


2.5 


21 


Go-3D.Uz-U.oo 


NH3(5,5) 


OIT^, 1 00 

Z.17 X l.zo 


2.6 


20 


r^QQft no n cq& 
LjOvSD.UZ-U.oo 


cont. 


VO \^ 1 QO 

Z. / Z X l.oZ 


n Q c; 
U.vSO 




-i cr no n 00 

Gc)45.UU-U.i2 


JN 113(4,4) 


1.14 X U.OO 


1.9 


43 


^ CT nn n 00 

Go45.UU-U. zz 


NH3(5,5) 


l.lZ X U.o4 


2.0 


37 


Ljo4o.UU-U.ZZ 


cont. 


n vft \y n Q/1 
U. /D X U.o4 


1 ft 
l.D 




Gool.77-U.o4 


M TU (A A\ 

JN 113(4,4) 


1.11 X 0.57 


4.0 


78 


/^Qt^i '7'7 n c; /I & 
LjOOI. / /-U.04 


1\.T1U (A A\ 
INII3 (4,4 ) 


ftft \y 1 01 

z.bD X l.zl 


2.8 


120 


LjoOI. ( /-U.04 


IN II3 (0,0) 


1 nn \y n c^fi 
l.Uy X U.OD 


A 1 
4. 1 


ftQ 


LjoOI. ( (-U.04 


1N113 (0,0) 


c;n \y 1 1 n 

z.oy X i.iy 


1 
Z. 1 


1 on 
IZU 


LjoOI. ( (-U.04 


cont 


^ ^ ^ \y c\ ftQ 
1.11 X U.Do 


1 1 
1.1 




G351.77-0.54 


cont 


000 , "1 r'r' 

vi.08 X 1.55 


1.8 




G0.55-0.85'= 


NH3(4,4) 


2.10 X 0.75 


2.6 


44 


00.55-0.85'= 


NH3(5,5) 


1.81 X 0.54 


2.5 


32 


00.55-0.85'= 


cont 


1.46 X 0.54 


1.3 


23 


019.47-0.17'= 


NH3(4,4) 


4.51 X 0.85 


3.1 


17'' 


019.47-0.17'= 


NH3(5,5) 


6.48 X 1.18 


4.8 


17'' 


019.47-0.17'= 


cont 


11.92 X 1.05 


0.35 


6 


118151-1208 


NH3(4,4) 


2.60 X 0.58 


2.9 




118151-1208 


NH3(5,5) 


2.53 X 0.58 


4.1 




118151-1208 


cont 


2.60 X 0.58 


0.24 


1.22 



For the line rms we used channel separations of 0.8 km s ^ . 

* Only antennas 1 to 5 were used for these lower-resolution images. 

'= Only limited shorter baseline ranges were used to produce these 

images. 

'' Peak flux of integrated image from 17 to 23kms~^. 
*= Negative because in absorption. 

This so urce was part of the Class II CH3OII maser 
sample by I Caswell et al.l (|1995f) . It is one of the strongest 
NH3(4,4) and (5,5) emiting sources in our sample (Table 
[2]) . With this strong emission it is also easy to detect and 
map the hyperfine structure components. Figure [7] shows 
the 1st moment maps of the main hyperfine components 
of the NIl3(4,4) and (5,5) lines excluding antenna 6 to 
better show the large-scale velocity gradient. We clearly 
identify a velocity gradient in north-west south-east di- 
rection centered around the CH3OII maser position. 

Including now antenna 6 in the imaging process to also 
study smaller-scale sub-structure, Figure [5] presents the 
1st moment maps of the main hyperfine components of 
the NH3(4,4) and (5,5) lines as well as the 1st moment 
map corresponding to the most blue-shifted hyperfine 
structure line (offset by 2.45 MHz from the main line). 



While the most blue-shifted hyperfine line exhibits a rela- 
tively clear velocity gradient in approximately northwest- 
southeast direction, the picture is less clear for the main 
hyperfine structure line of both NII3 transitions. Inspec- 
tion of the spectra shows that the hyperfine satellite lines 
have almost the same intensities as the central main line 
which indicates the extremely high optical depth. With 
such high optical depth we only see the outer layers of 
the rotating and likely collapsing core. In these lines 
we also identify a velocity gradient, however, the highest 
velocities are not found at the southeastern edge as for 
the blue-shifted NIl3(4,4) hyperfine line, but it is shifted 
about 1" inward of the core. This feature is reminiscent 
to the so-called bull 's-eye structu r e obse rved in NIl3(3,3) 
absorption lines bv lSollins et al.l (|2005 ^ toward the very 
luminous high-mass star-forming region GlO.6-0.4. The 
bull's-eye structure features the g as with the highest red- 
shift with respect to the wisr, and lSolhns et al.l (|2005f ) in- 
terpret this structure as caused by spherical infall of the 
core. Similarly, in our source G327.3-0.60, the most red- 
shifted feature at -42kms~"'^ with respect to the uisr of 
about — 46kms^^ (Table [1] and Figure [8] left and middle 
panel) exhibits the similar structure where the surround- 
ing gas in all directions features lower velocities again. 
Therefore, we may witness here some more spherical in- 
fall motions in the outer core that is traced by the more 
optically thick main hyperfine components as well. To 
further emphasize the complex morphological and kine- 
matic structure of this core. Figure [9] presents a chan- 
nel map of the region. We clearly see several distinct 
clumps which may all be associated with infall motions. 
At the near kinematic distance of ^3.1 kpc (Table [1]), 
the spatial extent of the structure (~5.5") corresponds 
to an approximate linear extent of 17000 AU. A position- 
velocity cut of the NIl3(4,4) main hyperfine line through 
the bull's-eye feature with a position angle of 140 de- 
grees east of north (Figure [10] left panel) also highlights 
the complexity of the structure without any Keplerian 
disk-like signature. 

The picture appears slightly different if one investigates 
the more optically thin blue-shifted NH3(4,4) hyperfine 
line (Figure [8] right panel). There we do not see the 
bull's-eye structure, but the 1st moment map features a 
more consistent velocity gradient in northwest-southeast 
direction. Since we do not know for certain about any 
outflow structure, assigning this velocity gradient to a 
rotating disk-like entity is dubious. Nevertheless, since 
these highly excited lines are not expected to trace out- 
flow motions, it is tempting to interpret the NIl3(4,4) 
hyperfine line 1st moment map as tentative evidence for 
rotation in this core. 

Combining the previous spherical infall signature from 
the more optically thick main hyperfine structure lines 
with the tentative rotational signature of the more op- 
tically thin blue-shifted hyperfine structure line, these 
data are consistent with a large-scale infalling envelope 
that appears spherical in the outer regions and exhibits 
stronger rotation signatures further inside due to the con- 
servation of angular momentum. 

Figure[8]also presents the corresponding line width 2nd 
moment maps. In all three cases the line width distri- 
bution peaks approximately in the center of the rotating 
structure, offset from the previously discussed bull's eye 
feature, and closer to the position of the Class II CII3OII 
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maser features. This again indicates that the hkely cen- 
ter of star formation activity, and hence the position of 
the central protostar, is actuaUy at the center of the core 
and not toward the bull's eye position. 

We also detect a 1.25 cm continuum source within our 
field, however, this is again more than 10" offset from our 
target and can therefore be considered in our context as 
unrelated. 

4.1.5. G328.81+0.63 (IRAS 15520-5234) 

As shown in Figure [T^ this region exhibits a small 
cluster of 1.25 cm continuum sources with more than 10 
emission peaks within the inner 21000 AU. The Class II 
CII3OII masers are located at the edge of the UCHii re- 
gion. The 1.25 cm peak flux is 136 mJy beam" ^. While 
we do not detect NH3 emission from the central re- 
gion, both NH3 transitions show strong absorption fea- 
tures toward the strongest cm continuum peaks. Fig- 
ure [12] shows the 1st moment maps of these absorption 
features. While the easternmost peak exhibits the ab- 
sorption peak around -39.3 km s~^, going further to the 
east we find peak velocities around -45.9, -44.8 and - 
44.0 km s""'^. Hence, while there are clear velocity differ- 
ences between the sources, there is no consistent velocity 
gradient across several sub-peaks. A closer inspection of 
the velocity structure of individual sub-peak reveals that 
- si milar to sources G3 27.3-0.60 (Section |4T11) or G10.6- 
0.4 ([Sollins et al.ll2005[ ) - toward the strongest eastern cm 
continuum source the absorbing gas toward the peak is 
systematically red-shifted compared to the edge of the 
peak position (Figure [T5)) . In the above discussed sce- 
nario fSection I4.1.4p . this again indicates approximately 
spherically infalling gas toward the Hii region. Follow- 
up observations of hydrogen recombination lines would 
be required to test whether the gas continues to infall 
through the Hii region, and the source may hence still 
be accreting, or whether it is stopped before an d the 
accre tion processes have already terminated (e.g., iKetol 
[2OOI) . 

The line width distributions of the absorption features 
are shown in Figure 1141 and we clearly see a line width 
increase toward the continuum peak positions. There are 
additional line width maxima in the Figure between the 
two main continuum peaks and at the northern edge of 
the strongest feature, corresponding to the most negative 
velocity features seen in Figure [T51 However, we refrain 
from a further interpretation of these edge effects. 

It is interesting to note that although we do not de- 
tect NH3 in emission from the region covered by the 
UCHii regions, we do indeed detect strong NH3 emission 
in both transitions on larger scales around the UCHii 
region. Figure fTSl shows the integrated NH3(4,4) emis- 
sion of the region at lower spatial resolution (we only 
used the baselines between and 45 fcAwith a synthesized 
beam of 6.1" x 4.1"), and we clearly identify a ring-like 
structure with a diameter exceeding 20". Since these 
are interferometer observations filtering out the largest 
spatial scales, it is likely that the real emission is even 
more extended. While extended NH3 emission around 
forming high-mass stars has previously been observed, 
for example toward NGC 6334 1(N) it was only seen in the 
low-energy (1,1) and (2,2) transitions, whereas the high- 
energy transitions up to the (6,6) line were observed onl y 
toward the central sources (|Beuther et al.l[2005l[2007bl) . 



Therefore, it is surprising that we witness large-scale hot 
NH3 emission around G328.8H-0.63 from lines with ex- 
citation temperatures as high as 295 K (see section [T]) . 
This implies that the central formed/forming high-mass 
stars have heated up significant amounts of gas out to 
distances exceeding 30000 AU from the center to tem- 
peratures in excess of 100 K without yet destroying the 
surrounding gas envelope. This may be interpret ed as 
support for the proposal bv lKrumholz et aTl (|2007b[ ) that 
radiative feedback from the central protostar is able to 
heat up the surrounding envelope strong enough that fur- 
ther thermal fragmentation will be largely suppressed. 

Figure [T5I also presents a position- velocity cut through 
the center of the large-scale map in east-west direc- 
tion. While the main emission is approximately be- 
tween -40 and -42kms^^, one absorption peak is around 
-40kms^^ while the rest of the absorption is more blue- 
shifted < — 44kms~^. Correlating these features with 
the higher-resolution absorption maps in Figure [121 the 
-40kms~"'^ component belongs to the eastern absorption 
peak, the features around -44 to -45kms~^ to the rel- 
atively extended east-west continuum ridge and the ab- 
sorption feature at ~ — 46kms~^ to the 2nd strongest 
absorption weak (2nd peak from west in Fig. [T^ . While 
for the strongest absorption peak relative motions with 
respect to the ambient cloud are not distinguishable, the 
blue-shifted absorption data for the 2nd strongest ab- 
sorption feature are indicative for expansion motion of 
the molecular envelope. 

4.1.6. G331. 28-0.19 (IRAS 16076-5134) 

The 1.25 cm continuum maps reveal 2 sources, how- 
ever, the NH3 emission is only associated with the weaker 
northern source (Fig. [T6)) which is also the Class II 
CH3OH maser emitter. The various outfiow tracers like 
H2 emission, Spitzer "green fuzzies" , SiO emission or 
7 mm continuum emission all indicate an outflow direc- 
tion approximately in northeast-southwest direction (Ta- 
ble [T]). In contrast to that, our NH3(4,4) and (5,5) ob- 
servations are more indicative of a velocity gradient in 
northwest-southeast direction, approximately perpendic- 
ular to the outfiow. It is interesting to note that the Class 
II CH3OH maser orientation appears to align closely with 
the NH3 emission. One peculiarity of the NH3 data is 
that the velocity gradients of the (4,4) and (5,5) lines 
are in approximate opposite directions. Since we do not 
identify clear Keplerian-like signatures, we do not ob- 
serve a real accretion disk, however, the orientation of 
the NH3 velocity gradient is strongly suggestive of ro- 
tating material perpendicular to the outflow axis. This 
gas may feed an accretion disk closer to the center of the 
core. 

Figure [16] also shows the 2nd moment line width dis- 
tribution. While for the (4,4) fine the data are less clear, 
a broader line width toward the center close to the Class 
II CH3OH maser features is observed for the (5,5) line, 
consistent with a central location of the protostar and 
hence the center of active infall. 

4.1.7. G336.02-0.83 (IRAS 16313-^840) 

Except for a potential Class II CH3OH maser ve- 
locity gradient in approximately north-south direction 
() Walsh et al.l Il998f ) , little else is known about this re- 
gion. We do detect weak NH3(4,4) and (5,5) emission 
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associated with the Class II CII3OII masers, however, 
only weakly when excluding the long baselines associated 
with antenna 6 from the data reduction. Figure [T7] shows 
the 1st moment maps. Although the spatial resolution 
does not allow to identify an obvious velocity gradient, in 
particular the NIl3(5,5) 1st moment map is indicative of 
a potential velocity gradient in the north south direction, 
parallel to the Class II CH3OH maser features. The 2nd 
moment line width distributions also shown in Figure [T71 
does not exhibit a prominent line width signature, pre- 
venting us from further interpretation. Furthermore, we 
do detect 1.25 cm continuum emission from the region, 
however, the peak is approximately 30" shifted to the 
north and can be hence considered as unrelated to the 
NII3 emission (Fig. fTS)) . 

4.1.8. G345.00-0.22 (IRAS 17016-4124) 

Class II CH3OH maser emission is det ected toward tw o 
positions approximately 4" apart fW alsh et al.l I1998D . 
While the eastern maser position is associated with 
1.25cm continuum emission and NH3(4,4) and (5,5) ab- 
sorption, the western maser peak is associated with 
NH3(4,4) and (5,5) in emission (Fig. [TH]). Although 
both maser groups appear to be spatially approximately 
aligned with an east-west orientation, they do not show 
a clear velocity gradient. This is at least different for 
the NH3 emission toward the western peak position 
which exhibits a clear velocity gradient in approximately 
east-west direction. Toward the eastern peak, most of 
the NH3 absorption is blue-shifted with respect to the 
visT ~ —26.8 (Table [1]), indicative of expanding gas. The 
2nd moment line width distribution (Figure I20p toward 
the western NH3 emission peak shows a line width in- 
crease approximately toward the central Class II CH3OII 
emission features. The signatures toward the eastern ab- 
sorption features are less conclusive. 

Therefore, while the NH3 absorption and 1.25 cm con- 
tinuum emission toward the eastern peak is consistent 
with an expanding UCHil region, the NH3 emission data 
toward the western peak position are indicative of a ro- 
tating structure that is likely still associated with ongo- 
ing high-mass star formation. 

4.1.9. G351. 77-0.54 (IRAS 17233-3606) 

This region is one of the previously most studied 
sources in our sample. It exhibits linear CII3OII 
maser features approximately in northeast-southwest di- 
recti on aligned with a CO outflow of similar orienta- 
tion (iNorris et al.lll993l: IWalsh et al.]|1998t iLeurini et all 
l2008f) . Furthermore, it exhibits an OH maser veloc- 
ity gradient in approximately north-south direction as 
well as a II2O maser structure that can be either inter- 
preted as a ring potentially associated with several re- 
cently identified cm sub-sources fFigs. [2T]and [22|) . or as 
well as a velocity g radient in northeast-southwest direc- 
tion (j^rsteretSllIllS; iFish et al.l 120051 : IZapata et all 
|2008[) . We do detect strong NH3(4,4) and (5,5) emis- 
sion from the CII3OH maser position with a clear veloc- 
ity gradient in ESE-WNW direction, approximately per- 
pendicular to the outflow and maser orientation. This 
signature can be depicted in the lower-resolution image 
excluding antenna 6 in the data reduction to highlight 
the larger-scale rotating signature fFig. [^Tjl . as well as in 
the highest-resolution images including antenna 6 to also 



show smaller-scale sub-structures (Figs. [22l fc [23|. This 
can be interpreted as good evidence of rotational motion 
of the core. Although the channel map in Fig. [221 exhibits 
a clumpy structure of the rotating gas similar to G327.3- 
0.60 (section [4. 1.4p . the velocity gradient can also clearly 
be depicted in the clumpy substructure. While this may 
not be significant, it is intere sting to note t hat th e VLA 
cm continuum sources from IZapata et al.l ()2008f ) corre- 
late in some spectral channels with the clumpy molecular 
sub-structure (Fig. The position- velocity diagram 

in Figure [24] also shows this velocity gradient, however, 
again the structure is very clumpy and does not resem- 
ble what one would expect from a Keplerian disk. It 
more resembles a large-scale rotating envelope structure 
that may feed the potential inner accretion disk. This 
is consistent with the 2nd moment line width distribu- 
tion (Fig. [22]) which is centrally peaked and hence consis- 
tent with increasing rotational velocities toward the cen- 
ter. This large-scale rotating envelope has a projected 
diameter of ^ 5" corresponding at the given distance 
of '^2.2kpc to a an approximate extent of ^^llOOOAU. 
Since this structure encompasses all 6 cm continuum 
sources identified bv IZapata et al.l (|2008D . this rotating 
envelope may even feed several smaller independent ac- 
cretion disks that could be associated with individual 
sub-sources. Nevertheless, since the outflow and rotation 
signatures do not appear very disturbed by the multiplic- 
ity, it is interesting that the general outflow and perpen- 
dicular rotation structures appear to be dominated by 
one object. 

We also detect strong 1.25 cm continuum emission from 
a nearby ultracompact Hii region separated by about 12" 
to the east (Fig. [^5)1 . At the highest spatial resolution 
available, this UCHii region splits up into 4 sub-sources. 
In addition to the UCHii region, we detect weak cm con- 
tinuum emission with a peak flux of ~14mJybeam~^ 
toward the NII3 emission source. This cm structure is 
not compact but rather extended and hence neither re- 
sembles a jet-like feature nor a hypercompact Hll region. 

4.1.10. GO. 55-0.85 (IRAS 17470-2853) 

This is again one of the sources with little additional in- 
formation. The region hosts two Class II CH3OII maser 
sites separated by ~ 3" but no obvious velocit y gradi- 
ents are present within them (jWalsh et al.lll998f) . We do 
detect the NH3(4,4) and (5,5) emission from both maser 
positions, and Figure [26] presents the corresponding 1st 
moment maps of the spectra lines. The two Class II 
CII3OII maser feature are clearly connected in the NH3 
emission but neither the (4,4) nor the (5,5) transition 
exhibits any conspicuous velocity structure. However, 
the 2nd moment line width distribution (Fig. 1261 shows a 
double-peaked structure, where the two line width peaks 
are associated with the two Class II CH3OH maser fea- 
tures and the peaks of integrated NH3 emission (Fig. Wf\ . 
This is further evidence for active star formation activity 
associated with both emission peaks. 

Furthermore, we identify 1.25 cm continuum emission 
from a likely associated UCHil region directly south of 
the NH3 thermal and CH3OH maser emission (Fig. Wf\ . 
In addition to this, we identify a second 1.25 cm con- 
tinuum peak at the ^ 6a level clearly associated with 
the NH3 peak emission and the northern CH3OH maser 
position. 
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4.1.11. GI9.47+O.I7 (IRAS 18232-1154) 

This region exhibi ts a CO outflow in a. pproximate 
NNE-SSW direction (jLongmore et al.|[2007al Longmorc 
et al. in prep.). We do detect both NH3 hnes close to 
the Class II CH3OH maser position of this region as 
well (Tabled]), however, the signal-to-noise ratio is rela- 
tively poor, prohibiting good quality moment maps and 
hence deriving reliable velocity gradients. Furthermore, 
we do detect extended 1.25 cm continuum emission to- 
ward the north of the NH3 emission peak. Figure [551 
presents an overlay of the 1.25 cm continuum emission 
with the integrated NIl3(4,4) emission. The integrated 
NIl3(4,4) emission map shows several additional Aa fea- 
tures distributed in the vicinity which may indicate more 
extended NH3 emission just f iltered out by our inter fer- 
ometer observations (see also iLongmore et "al1l2007al l. 

4.1.12. IRAS 18151-1208 

As outlined in section [3l this region does not satisfy 
the selection criteria of the rest of the sample. However, 
since it is also a Class II CH3OH maser source and ex- 
hibits str ong outflow and d isk signatures (|Beuther et al.l 
I2002cl lbl; iDavis et all |2004 Fallscheer et al. in prep.), 
we considered it a good addition for this observing run. 
However, unfortunately it remained undetected in both 
NH3 transitions. Nevertheless, we did detect at the 5a 
confidence level 1.25 cm continuum emission associated 
with the CH3OH Class II maser peak (Fig.[29]and Table 

4.2. Spectral fitting and temperature determination 

The low transition NH3 lines (e.g., the (1,1) and (2,2) 
lines) are known to be an excellent thermometer for the 
cold components of the gas wit hin molecular clouds (e.g., 
IWalmslev fc UngerechtsI [l983[ ). Similarly, we may be 
able to use the high-transition lines here to estimate the 
temperatures of the warm gas observed in these regions. 
The advantage of NH3 is that one observes the whole 
hyperfine structure simultaneously and hence should be 
able to derive the optical depth of the lines. In the case of 
the NH3(4,4) and (5,5) transitions, the relative intensi- 
ties in the optically thin regime of the satellite lines with 
respect to the main central hyperfine components are ap- 
proximately 2 and 1%, respectively. Figures \30\ [311 and 
[32] show example spectra of each source and our attempts 
to fit the whole hyperfine structure to derive the spectral 
parameters. However, in most cases the fits do not rep- 
resent the data well. This is mainly due to the extremely 
high optical depth where the satellite lines reach about 
the same intensities as the main central component. In 
such a case the fits give flat-topped spectra which are 
not observed. This discrepancy indicates that there is a 
temperature gradient along the line of sight which is not 
taken into account by the fitting procedure. More ad- 
vanced radiative transfer calculations would be required 
to reproduce the spectral shape which is out of the scope 
of this paper. Therefore, we are not able to get accu- 
rate temperature estimates for the target sources. How- 
ever, based on the high excitation temperatures of the 
two lines (i?iowcr of 200 and 295 K, respectively) and the 
high observed brightness temperatures between several 
10 and more than 100 K in most sources (see Figs. [3011311 
and [5^ , it is reasonable to assume that the average gas 
temperatures in the observed regions exceeds 100 K. 



While temperature gradients increasing toward the 
central protostars do also increase the thermal line width 
toward the center, this is unlikely to explain the central 
line widths increases as observed in several of the 2nd 
moment maps. The thermal line width scales with the 
square-root of the temperature (Authcim cx: T*^'^). As- 
suming a temperature gradient T{r) oc r"^ '*, the ther- 
mal line width scales like Awthcrm oc r""'^. For example, 
assuming 100 K temperature at a core edge with a radius 
of 5000 AU the thermal line with Authcrm is ~ 0.5km~^. 
With the above relation Awthorm at an inner radius of 
500 AU should be ~0.8kms-i. If we take G327.3-0.60 
in Fig. [8] as an example, the central values exceeding 
4kms~^ are far broader than any thermal line broad- 
ening could produce. The associated velocity gradients 
suggest that the dominating reason for the observed line 
width broadening should be due to rotation. 

4.3. Morphologies of extended continuum emission 

Three of the regions we have imaged show clear evi- 
dence for extended continuum emission - G305. 21-1-0. 21, 
G316.81-0.06 and G336.02-0.83. In aU cases, the con- 
tinuum emission appears to resemble a simple circular 
shape of diameter 30 — 40". We believe that in each 
case, extended continuum emission is present, however, 
we question the validity of the morphologies shown in 
the continuum images. This is because an angular size of 
^ 100" corresponds to a baseline length of approximately 
2.5 kA, which closely matches the shortest baseline used 
in these observations (31m). The second shortest base- 
line of 199 m corresponds to spatial scales of ~ 15". 
While scales below 15" are sampled relatively well by our 
observations, only a single baseline covers larger spatial 
scales, and structures above ~ 100" are completely fil- 
tered out. Therefore, our observations are sensitive to 
extended structures over only a small range of sizes and 
thus do not represent the true morphology of extended 
emission. 

Recent observations by ILongmore et al.l (|2009[ ) of 
G316.81-0.06 at 18.8 GHz show two sources on either side 
of the emission we show in Figure [5] These recent obser- 
vations include a better sampling of the UV plane over 
sca les corresponding to th e extent of the emission found 
by ILongmore et al.l ()2009l ) and so are a more accurate 
representation of the emission in this region. To double- 
check, by excludi ng baselines > 4 kA from the continuum 
data presented in ILongmore et aTl (|2009| ). we recover an 
image similar to that in the current paper. 

We therefore conclude that the morphology of con- 
tinuum emission in G305.21-h0.21, G316.81-0.06 and 
G336.02-0.83 is probably not accurately represented by 
that shown in Figures [3] [5] and [18] respectively. We cau- 
tion the reader on interpretation of extended emission 
with similar interferometer configurations. 

5. GENERAL IMPLICATIONS 

Table [3] summarizes the general results regarding ro- 
tational signatures for the whole sample. Except for the 
source IRAS 18151-1208 which did not satisfy the origi- 
nal selection criteria from the rest of the sample, all other 
sources were clearly detected in the high-excitation NH3 
lines, implying that our sample selection criteria were 
well chosen. Out of the remaining eleven sources, six 
show signatures of rotation and/or infall which can be 
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considered as strong evidence of ongoing high-mass star 
formation activity. While the rotational signatures vary 
between typical velocity gradients perpendicular to the 
outflows to more spherical infall signatures and infall sig- 
natures from absorption lines, we do not find clear signs 
of Keplerian rotation. This implies that although we 
have achieved very high angular resolution, mostly bet- 
ter than 1" (Table [2]), the corresponding linear scales 
(Tables [T] & [2|) do not follow typical disk-like Keplerian 
signatures. Since outflows and jets are known for most of 
the sources, and sinc e these are believ ed to be accelerated 
via disk- winds (e.g.. lArce et al.l[2007[ ). the corresponding 
accretion disks should be smaller. Furthermore, the spec- 
tra presented in Figs. [501 [5T] and[5^show that even highly 
excited lines still have a very high optical depth. Hence 
we are only seeing the r = 1 surface of these spectral 
lines and do not trace the innermost regions. Hence the 
accretion disks can still be obscured by the high optical 
depths. 

Compared to other examples where s patial fla ttening is 
observed on scales > 10" AU (e.g., MlT lChini et al. 2004, 
IRDC18223-3 FaUscheer et al. 2009), it is surprising that 
we do not find a single source where spatial flattening is 
observed in our data. To the contrary, in several sources 
where we see clear velocity gradients indicative of rota- 
tion, the spatial structure of the gas is distributed over 
scales > 10"* AU in a very clumpy fashion (e.g., G327.3- 
0.60 or G351.77-0.54). While we find infaU signatures 
in absorption lines (e.g., G328. 81-1-0. 63), it is also inter- 
esting to note that in at least one source (G327.3-0.60) 
we observe the "bull's-eye" signature in the 1st moment 
map that was attr ibuted to spherical infall motions by 
ISolUns et all ([2001 for one of the most high-mass star- 
forming regions GlO.6-0.4. 

The high-excitation lines are always found toward the 
CH3OH Class II maser positions, reinforcing the idea 
that these masers are excellent tracers for early stages 
of high-mass star formation. However, the data do not 
allow us to draw clear conclusions whether these maser 
are disk or outflow associated. 

In addition to the NH3 data, toward most fields we de- 
tect 1.25 cm continuum emission from either UCHii re- 
gions associated with our targets, or, more typically, ex- 
tended continuum emission offset more than 10" tracing 
more evolved Hil regions in the field of view, which how- 
ever is not directly associated with our target sources. 
This highlights that high-mass star formation rarely pro- 
ceeds in isolation but that other potentially more evolved 
regions often can exist within the same star formation 
complex.. 

6. CONCLUSION 

We present a large observing campaign to search for 
rotational motions in high-mass star-forming regions 
via high-excitation NH3 line observation with high 
spatial resolution. While we detect all sources from 
the original sample (excluding the one source that was 
chosen via other selection criteria) , in more than 50% of 
them rotational and/or infall motions can be identified. 
Several more general conclusions can be drawn from this 
observing campaign. 

(a) High-excitation NH3 lines are very good tracers of 
the dense gas within hot-core type young high-mass 



star- forming regions. In this sample, the NH3(4,4) and 
(5,5) emission is always associated with the CH3OH 
Class II methanol maser emission, which is a well-known 
signpost for high-mass star formation. 

(b) We identify rotational/infalling motions in half of 
the observed sources. 

(c) The signatures comprise simple velocity gradients as 
well as more spherical infall signatures. 

(d) Although the spatial resolution is of order a few 
1000 AU, we do not identify obvious Keplerian signa- 
tures. 

(e) While high optical depth is an issue, the data never- 
theless indicate that accretion disks resembling those o f 
low- mass star- forming regions (e.g., ISimon et al.l |2000| ) 
have to be typically of smaller sizes, not yet resolved 
by our observations. Likely such smaller accretion disks 
could be fed by the rotating envelopes observed here. 

(f) In many fields of view we do detect additional Hii 
regions that are not directly linked to our targets but 
that reside in the same projected area on the sky. This is 
additional evidence for the clustered mode of high-mass 
star formation. 

Where to go from here? Clearly there are differ- 
ent paths to be followed. While our approach for this 
study was biased toward hot-core type sources based 
on our spectral line selection criteria, other observa- 
tional approaches can target sources in different evo- 
lutionary stages as well as of different luminosities. In 
particular, observations at mm wavelengths with broad 
spectral bandpasses allow us to observe several spectral 
lines simultaneously which can trace rotati onal motions 
in sources of different characteristics (e.g . , iZhangj 120051 : 
iCesaroni et "all 120071: iBeuther et all 120091 ) . For example, 
Qizhou Zhang and collaborators started a similarly large 
project searching for disk signatures toward a sample of 
young high-mass star-forming regions with the Submil- 
limeter Array, and they also detected rotational signa- 
tures toward all targets (Qizhou Zhang et al., in prep.). 
While as outlined in section [5] flattened structures do ex- 
ist on scales exceeding 10* AU, most data indicate that 
the Keplerian accretion disks should reside on smaller 
scales, probably < 1000 AU. This is also shown in the re- 
ce nt 3-dimensional radiative hy drodynamic simulations 
bv lKrumholz et all ()2007bl . [2009l ) where most of the disk- 
like structures are contained with a radius of ~500AU. 
Furthermore, recent work has also shown that accretion 
disks exceeding 150 AU are prone to fragmen tation (e.g., 
iKratter k Matznerj|200l I Vaidva et al.ll2009f) . hence, it is 
not even expected to find Keplerian signatures on much 
larger scales. 

While we can already achieve tremendous progress 
with the existing instrumentation, clearly we are still 
lacking spatial resolution if we want to understand the 
structure and physical processes of the more proper cen- 
tral accretion disks. Future instruments like ALMA and 
eVLA will allow us to make considerable progress in this 
direction for the cold gas components. Furthermore, 
due to the large accretion rates, strong viscous forces 
and central luminous sources, a considerable fraction of 
warm/hot gas is expected to exist with the high-mass 
disks. Therefore, imaging these warm/hot dust and gas 
components at high spatial resolution with future instru- 
ments like MIRI on JWST will allow us to constrain the 
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TABLE 3 

Summary results for rotational signatures. 



Source 



Results 



G305.21+0.21 
G316. 81-0.06 
G323.74-0.26 
G327.3-0.60 
G328.81+0.63 
G331.28-0.19 

G336.02-0.83 
G345. 00-0.22 
G351. 77-0.54 
GO.55-0.85 
G19.47-0.17 

properties of thes e MerspaBsaits iiy<}»]ai6b;<aieflLe detail. It 
is hkely that only a concerted effort at long and short 



Rotating structure perpendicular to outflow, no Keplerian signature. 
Inconclusive. 
Inconclusive. 

Infall and rotational signatures. Clumpy large structure. 
Infall toward UCHii regions. Large-scale hot gas in emission. 
Velocity gradient perpendicular to outflow suggestive of rotation and infall. 
However, differences between lines. 
Inconclusive. 

Rotation in western peak, expanding UCHll region in eastern peak. 
Rotating structure perpendicular to outflow, no Keplerian signature. 
Inconclusive. 
Inconclusive. 

transfer modeling will give us a real understanding of the 
working of high-mass accretion disks. 



wavelengths combined with hydrodynamic and radiation 
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Fig. 1. — G305.21+0.21 intensity weighted velocity (1st moment, top row) and line width (2nd moment, bottom row) maps of the 
main hyperfine com ponents of NH3(4,4) (left) and NH3(5,5) (right). The triangle marks the main Class II CH3OH maser position from 
INorris et all l(T993h . the arrows in the top-left panel outline the approximate direction from the assumed outflow (e.g., De Buizer 2003i). 
the synthesized beams are shown at the bottom-left of each panel, a scale-bar is presented in the top-left panel, and the 0/0 position is 
given in Table [Tl 
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Fig. 3. — 1.25cm continuum emission in the field of G305. 21-1-0. 21. The contour levels are in 3cr steps (Tabled. The star marks the 
position of our primary NH3 target within the field, and a scale-bar is shown at the bottom-right. 



Fig. 4. — G316.81-0.06 intensity weighted velocity (1st moment, top row) and line width (2nd moment, bottom row ) maps of the mai n 
hyperfine components of NH3(4,4) (left) and NH3(5,5) (right). The triangles mark the Class II CH3OH maser positions ll Walsh et aLllTOOSl ). 
and the synthesized beams are shown at the bottom-left of each panel, a scale-bar is presented in the top-left panel, and the 0/0 position 
is given in Table [T] 
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Fig. 5. — 1.25cm continuum emission in the field of G316.81-0.06. The contours start at the 5cr level and continue in 3cr steps {Table[2]l. 
The star marks the position of our primary NH3 target within the field, and a scale-bar is shown at the bottom-right. 



Fig. 6. — G323.74-0.26 intensity weighted velocity {1st moment, top row) and line width {2nd moment, bottom row) maps of the main 
hyperfine components of NH3{4,4) {left) and NH3(5,5) (right). The triangles mark the Class II CH3OH maser positions (Walsh et al. 1998), 
the synthesized beams are shown at the bottom- left of each panel, a scale- bar is presented in the top-left panel, and the 0/0 position is 
given in Table [Tl 



Fig. 7. — G327.3-0.60 intensity weighted velocity {1st moment) maps of the main hyperfine components of NH3{4,4) (left) and NH3(5,5) 
(right) excluding antenna 6 to better show the larger-s cale structure . The triangles mark the positions of the Class II CH3OH masers 
(James Caswell, private communication as referred to bv lCaswelllll99^ . the synthesized beams are shown at the bottom-right of each panel, 
a scale-bar assuming near kinematic distance is presented in the left panel, and the 0/0 position is given in Table[l] 



Fig. 8. — G327.3-0.60 intensity weighted velocity {1st moment, top row) and line width (2nd moment, bottom row) maps of the main 
hyperfine components of NH3{4,4) (left) and NH3{5,5) (middle) including also antenna 6 data to also show the smaller-scale sub-structure. 
The right panel shows the 1st moment map of the most blue-shifted NH3(4,4) hyper fine c omponent. The triangles mark the positions of 
the Class II CH3OH masers (James Caswell, private communication as referred to bv lCaswel]|il998i) . the synthesized beams are shown at 
the bottom-left of each panel, a scale-bar assuming near kinematic distance is presented in the bottom- left panel, and the 0/0 position is 
given in Table [Tl 
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Channel map of the main hyperfine component of NH3(4,4) with a spectral resolution of 2kms^^ in G327.3-0.60. The contour 



Fig. 9 

levels (positive full lines, negative dashed lines) are in So" steps with a Icr value of 1.6mjybeam~^ 
the bottom-left of each panel, and the 0/0 position is given in Table[T] 



The synthesized beams are shown at 




-45 -40 
Velocity [km s ^] 

Fig. 10. — G327.3-0.60 position velocity diagram of the main hyperfine components of NH3(4,4) (left). The diagram is centered at offset 
(—12.5"/ — 9.1") with a position angle of 140 degrees from north (northwest-southeast direction). 
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Fig. 11. — 1.25 cm continuum emission in the field of G327.3-0.60. The contour levels start at 3a" and continue in 2a steps (Table|2l full 
lines positive, and dashed lines negative features). The star marks the position of our primary NH3 target within the field, and a scale-bar 
is shown at the bottom-right. 



Fig. 12. — G328.81+0.63: The grey-scales present the 1st moment maps of the absorption features of the main NH3(4,4) and (5,5) 
hyperfine components (top and bottom panels, respectively) against the ultracompact Hll region outlined in contours from the 1.25 cm 
continuum emissio n. The contouring is done in 6fT steps of 16.8 mjy beam~^ . The triangles mark the Class II CH3OH maser positions 
l|Walsh et al.|[T99^ . and the synthesized beams and scale-bars are plotted at the bottom left of each panel. 



Fig. 13. — G328.81+0.63 centered on the brightest eastern peak from Figure [T2l The grey-scale presents the 1st moment maps of the 
absorption features of the main NH3(4,4) hyperfine component against the ultracompact Hll region outlined in contours from the 1.25 cm 
continuum emission. The contouring is done in %a steps of 16.8 mjy beam~^ . 



Fig. 14. — G328. 81-1-0. 63: The grey-scales present the 2nd moment maps of the absorption features of the main NH3(4,4) hyperfine 
component against the ultracompact Hll region outlined in contours from the 1.25 cm continuum e mission. The contouring is done in 6(T 
stops of 16.8 mjy beam" ^ . The triangles mark the Class II CH3OH maser positions llWalsh et al.| [T998^. and the synthesized beam and 
scale-bar are plotted at the bottom left of each panel. 
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Fig. 15. — G328.81+0.63: Top: Integrated NH3(4,4) emission from -47.5 to -37.5 kms~^. To produce this map, we only used the short 
baselines between and 45 fcA. The contouring is done in itScr steps of 2.7mjy beam~^. The triangles mark the Class II CH3OH maser 
positions, the synthesized beam and the scale-bar arc plotted at the left. Bottom: The corresponding position-velocity diagram through 
the center in east- west direction. The contouring is done in 3(t levels of 6.6 mjy beam, measured in an 0.8 km si channel. Pull lines show 
emission, dotted lines absorption. 
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Fig. 16. — G331.28-0.19 intensity weighted velocity (1st moment, top row) and line width (bottom row) maps of the main hyperfine 
components of NH3(4,4) (left) and NH3(5,5) (right). The contours show the 1.25cm continuum emission with contour levels in 3a steps 
of 0.9 mjy beam~^ . The white triangles mark the Class II CH3OH maser positions ijWalsh et al.|[T993 ). the arrows in the top-left panel 
outline the approximate direction of the outflow (Table [TJ, the synthesized beams are shown at the bottom- left of each panel, a scale-bar 
is presented in the top-left panel, and the 0/0 position is given in Table[T] 



Fig. 17. — G336.02-0.83 intensity weighted velocity (1st moment, top row) and line width (2nd moment, bottom row) maps of the main 
hyperfine components of NH3(4,4) and NH3(5,5) excluding the long baselines associated with antenna 6. The triangles mark the Class II 
CH3OH maser positions by Walsh et al. (1998), the synthesized beams are shown at the bottom-left of each panel, a scale-bar adopting 
the near kinematic distance is presented in the top-left panel, and the 0/0 position is given in Table [l] 
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Fig. 18. — 1.25 cm continuum emission in the field of G336.02-0.83. The contour levels are at 3a- intervals (Table[2]l. The star marks the 
position of our primary NH3 target within the field, and a scale-bar adopting the near kinematic distance is presented in the bottom-right. 



Fig. 19. — G345.00-0.22 intensity weighted velocity maps (1st moment) of the main hyperfine components of NH3(4,4) and (5,5) in the 
left and right panels, respectively. The contours show the 1.25 cm continuum emission starting at the 4cr level and continuing in 6(t steps. 
While the main panels show the NH3 in e missio n, the inset presents the absorption against the continuum. The triangles mark the Class II 
CH3OH maser positions by Wa lsh et al.l 111998) . the synthesized beams are shown at the bottom-left of the two bottom panels, a scale-bar 
adopting the near kinematic distance is presented in the top-left panel, and the 0/0 position is given in Table [T] 



Fig. 20. — G345.00-0.22 intensity line width maps (2nd moment) of the main hyperfine components of NH3(4,4) and (5,5) in the left 
and right panels, respectively. The contours show the 1.25 cm continuum e mission starting at the 4a level and continuing in 6a steps. The 
triangles mark the Class II CH3OH maser positions bv lWalsh et al.l l|1998). the synthesized beams are shown at the bottom-left of the two 
bottom panels, a scale-bar adopting the near kinematic distance is presented in the top-left panel, and the 0/0 position is given in Table [Tl 



Fig. 21. — G351.77-0.54 intensity weighted velocity (1st moment) maps of the main hyperfine components of NH3(4,4) (left) and NH3(5,5) 
(right). The triangles mark the positions of the Class II CH3OH masers (Walsh et al. 1993), the arrows and star s in the left panel outline 
the approximate direction of the outflow (Table [TJ and the positions of the cm sources from lZapata et all II2008I) . The synthesized beams 
are shown at the bottom-left, a scale-bar is presented in the bottom-left panel, and the 0/0 position is given in Table [T] 



Fig. 22. — G351.77-0.54 intensity weighted velocity (1st moment, top row) and line width (2nd moment, bottom row) maps of the 
main hyperfine co mponents of NH3(4,4) (left) and NH3(5,5) (right). The triangles mark the positions of the Class II CH3OH masers 
llWalsh et al.| [l"993). the arrows and stars in the top-left panel outline the approximate direction of the outflow (Table [T} and the positions 
of the cm sources from Zapata ct al. (2008). The synthesized beams are shown at the bottom-loft of each panel, a scale-bar is presented in 
the top- left panel, and the 0/0 position is given in Table [T] 
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Fig. 23. — Channel map of the main hypcrfine component of NH3(4,4) with a spectral resolution of 2kms^^ in G351.77-0.54 . The 
contour levels (positive full lines, neg ative dashed lines) are in Scr steps with a la value of 3.3 mjy beam~^ . The crosses mark the cm 
continuum sources from lZapata et alj l l2008i) . and the synthesized beams are shown at the bottom-left of each panel, and the 0/0 position 
is given in Table [T] 
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Fig. 24. — G351.77-0.54 position velocity diagrams of the main hypcrfine components of NH3(4,4) (left) and NH3(5,5) (right). The 
diagrams are centered at Offsets (0"/0") with a position angle of 105 degrees from north (ESE-WNW). 
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Fig. 25. — 1.25 cm continuum emission in the field of G351.77-0.54. The top-panel shows the data excluding the long baselines associated 
with antenna 6, hence with lower spatial resolution (Table |2}. The contour levels are in 3cr steps (Table [2] full lines positive, and dashed 
lines negative features). The star marks the position of our primary NH3 target within the field, and a scale-bar is presented in the top 
panel. 



Fig. 26. — GO. 55-0. 85 intensity weighted velocity (1st moment, top row) and line width (2nd moment, b ottom row) m aps of the 
main hyperfine components of NH3(4,4) and (5,5). The triangles mark the Class II CH3OH maser positions by I Walsh et al.l (1998), the 
synthesized beams are shown at the bottom-right of each panel, a scale-bar adopting the near kinematic distance is presented in the top-left 
panel, and the 0/0 position is given in Tabic HI 
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Fig. 27. — The grey-scale with dotted contours shows the 1.25 cm continuum emission in the field of GO. 55-0. 85. The contour levels are at 
3<T intervals (Table [5J|. The full contours present the NH3(4,4) emission integrated from 13 to 21kms~^ with contour levels at 3cr steps of 
4.8mjybeam~^. The triangles mark the Class II CH3OH maser positions bv lWalsh et al.l I I1998I) . the synthesized beam of the continuum 
data is shown at the bottom left, and a scale-bar adopting the near kinematic distance is presented in the top-left. 
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Fig. 28. — The grey-scale with dotted contours shows the 1.25 cm continuum emission in the field of G19. 47-1-0. 17. The contour levels are 
at 3(T intervals (Table[2]l. The full contours present the NH3(4,4) emission integrated from 17 to 23kms~^ with contour levels starting at 
the 4(T level of 5.6mjybeam~^ and continuing at 3cr steps. The triangles mark the Class II CH3OH maser positions II Walsh et al-llTOOSl) . 
the synthesized beam of the NH3 data is shown at the bottom left, and a scale-bar is presented in the bottom-right. 
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Fig. 29. — 1.25cm continuum emission in the field of IRAS 18151-1208. The contour levels st art at the 3f7 level and continue in la- 
steps (Table[2J. The triangle marks the Class II CHaOH maser position from lBeuther et al.l Il2002d) . the synthesized beam is shown at the 
bottom left, and a scale-bar is presented in the bottom-right. 
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Fig. 30. — NH3(4,4) and NH3(5,5) spectra (left and right column) extracted toward the peak positions of the sources labeled in each 
panel. The histogram presents the data whereas the full lines show attempts to fit the whole hyperfine structure. Due to the very high 
optical depth, even this hyperfine structure fitting does not work well. 
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Fig. 31. — NH3(4,4) and NH3{5,5) spectra (left and right column) extracted toward the peak positio ns o f the sources labeled in each 
panel. For G328.81+0.63, the shown spectrum is extracted toward the large-scale western peak in Figure [TSl The histogram presents the 
data whereas the full lines show attempts to fit the whole hyperfine structure. Due to the very high optical depth, even this hyperfine 
structure fitting does not work well. 
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Fig. 32. — NH3(4,4) and NH3(5,5) spectra (left and right column) extracted toward the peak positions of the sources labeled in each 
panel. The histogram presents the data whereas the full lines show attempts to fit the whole hyperfine structure. Due to the very high 
optical depth, even this hyperfine structure fitting does not work well (it did not work at all for the (4,4) line of G351.77-0.54). 
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